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ON  THE  ZEROS  OF  EXPONENTIAL  POLYNOMIALS 

Cerino  E.  Avellar  and  Jack  K.  Hale 

ABSTRACT:  Suppose  r = , . . . , rM)  , r^  > 0,  Ykj  > 0 integers, 

k = 1,2,  ...,N,  j = 1,2,  Y^t  = Z.Ykjrr  The  purpose  of  this 

paper  is  to  study  the  behavior  of  the  zeros  of  the  function 

N " r 

h(A,r,a)  = 1 + Xj=1aje 

where  each  a^  is  a real  number.  More  specifically  if 

T(r,a)  = closure{Re  A:  h(A,r,a},  we  study  the  dependence  of  Z(r,a) 

on  r,a.  This  set  is  continuous  in  a but  generally  not  in  r. 

However,  it  is  continuous  in  r if  the  components  of  r are 

rationally  independent.  Specific  criterion  to  determine  when 

0 /£  TCr,a)  are  given.  Several  examples  illustrate  the  complicated 

nature  of  Z(r,a) . 

The  results  have  immediate  implication  to  the  theory  of 
stability  for  difference  equations 

M 

x(t)  - l A,  x(t-r,  ) = 0 
k=l  k K 

where  x is  an  n-vector,  since  the  characteristic  equation  has  the 
form  given  by  h(A,r,a).  The  results  give  information  about  the 
preservation  of  stability  with  respect  to  variations  in  the  delays. 
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The  results  also  are  fundamental  for  a discussion  of  the 
dependence  of  solutions  of  neutral  differential  difference  equations 
on  the  delays.  These  implications  will  appear  elsewhere. 
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ON  THE  ZEROS  OF  EXPONENTIAL  POLYNOMIALS 

by 

CERINO  E.  AVELLAR  AND  JACK  K.  HALE 

1.  Introduction.  Let  1R  = (-°°,°°),  1R*  = (0,«),  IR+  = 

a = (a^ , . . . , a^)  £IR  , r = (r^  , . . . , r^)  £ 0R  ) , Y j - (Yj^>..«,  Yj^)  » 

Y.v  nonnegative  integers,  j = 1,2, ...,N,  k = 1,2, 

JK  M. 

Y j • r = Yjkrk’  Our  purpose  in  this  paper  is  to 

of  the  real  parts  of  the  zeros  of  the  function 

N -Xy • *r 

(1.1)  h(X ,r,a)  = 1 + Ij=1a^e  J . 

More  specifically,  if 

(1.2)  Z(r,a)  = (Re  X:  h(X,r,a)  = 0} 

and  ZXr,a)  = cl  Z(r,a),  the  closure  of  Z(r,a),  we  study  the 
dependence  in  the  Hausdorff  metric  of  Z(r,a)  on  r,a.  It  is  shown  that 
7(r,a)  is  continuous  in  a with  a certain  type  of  uniformity  in  r. 

It  has  been  known  for  some  time  (see  Melvin  [6]  or  Henry  [4])  that 
Z\r,a)  is  not  continuous  in  r.  However,  we  show  that  it  is 
continuous  in  r if  the  components  of  r are  rationally  independent. 

We  also  give  a characterization  of  X(r,a)  in  a way  which  is 
amenable  to  computation.  For  the  case  in  which 
function  h(X,r,a)  is  given  as 


(1.3) 


2 


h(A,r,a)  = 1 + ^=1aje  •*  » 

the  characterization  of  Z(r,a)  is  more  complete  and  the  computation 
of  T(r,a)  can  be  given  rather  explicitly. 

Finally,  we  give  several  characterizations  of  the  property  that 
"Z(r,a)  n 1-6,6)  = <J> , 6 > 0;  that  is,  the  polynomial  h(A,r,a)  is 
hyperbolic.  The  case  "Z(r,a)  C C - <5)  , 6 > 0 is  also  discussed 
in  detail.  This  corresponds  to  uniform  asymptotic  stability. 

The  implications  of  the  results  for  difference  equations  are 
immediate.  In  fact,  consider  the  equation 

(1.4)  x(t)  - £jJ=lAkxtt"rk-)  = 0 

where  x e ]Rn  and  each  is  an  n * n matrix.  For  any 

<{>  e C = C([-h,  0]  ,]Rn)  , h >_  max{r^},  there  is  a unique  solution 
x = x(<J>)  of  (1.4)  for  t ^ -h  which  satisfies  x(4>)(t)  = <J>(t), 
t e t-h,0).  If  we  let  x(4»)(t  + 9)  = (S(t)4>)  (0)  , -h  < 6 < 0,  then 
SCt):  C -*■  C,  t > 0,  is  a strongly  continuous  semigroup  of  bounded 
linear  operators.  Furthermore,  if 

u(r,a)  = inf{b:  3 k with  |S(t)|  < ke'5t) 

then  it  is  known  (see  Henry  [4],  Hale  [2])  that 

oi(r,a)  = sup{Re  A:  h(A,r,a)  = 0} 
h(A ,r,a)  = det[I  - I^=1Akexp ( -Ark) ] . 


Therefore,  the  above  results  give  information  about  the  behavior  of 


the  order  ot(r,a)  of  the  semigroup  S(t)  as  a function  of  r,a. 

The  results  also  have  implications  for  neutral  functional 
differential  equations  of  the  type 

(1.5)  [x(t)  - lJJ=lAkx^t'rk-)1  = f(xt^ 

where  f:  C -*■  ]Rn  and  xt(®)  = x(t+®)>  -h  < 9 < 0.  The  solution 
operator  for  Equation  (1.5)  can  be  written  as  a sum  of  a completely 
continuous  operator  and  the  operator  S(t)  above  (see  Hale  [2]). 

If  f is  linear,  this  gives  information  about  the  spectrum  of  the 
solution  operator.  One  can  then  prove  certain  theorems  on  the 
continuous  dependence  in  the  delays.  Results  of  this  type  will  appear 
in  Avellar  and  Hale  [1] . 

2.  Continuous  dependence.  In  this  section,  we  present  some  results 
on  the  dependence  of  the  set  Z"(r,a)  on  r,a.  We  need  the  Hausdorff 
metric  which  is  defined  as  follows: 

For  any  sets  E,FC1R  and  any  point  p e IR,  let 

(i)  dCP.E)  - inf  | P-t I 
teE 

(ii)  fiCE.F)  = sup  d (P , F ) 
peE 

_ Ciii)  D (e , F ) = max{6(E,F),  6CF,E)}. 

The  number  D(E , F)  is  called  the  Hausdorff  distance 

l | 

between  the  sets  E,F  in  IR. 

We  need  the  following  result  fran  Levin  [14,  p.  268],  the 

proof  of  which  is  omitted. 
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Lemma  2.1.  For  a given  a < 0,  the  following  conclusions 
hold : 

(i)  There  is  an  integer  p such  that,  for  all  real  t, 
there  arc  no  more  than  p zeros  of  h in  the  box 

{X:  a < Re  X < 3;  t < ImX  < t + 1}. 

I 

(ii)  For  any  6 > 0,  there  is  an  mC<5)  > 0 such  that; 
whenever  a < Re  X < 3 and  X is  at  a distance  > 6 from 
every  zero  of  h,  one  has  |h(X)|  > m(,6). 

Our  first  objective  is  to  obtain  an  interval  which  con- 
tains Z(a,r).  Observe  that  X = u + iv  satisfies 
h(X,a,r)  =0  if  and  only  if 

n -uy.t  -ivVr  E , , -yVr  i(<*VvVr) 

0 = I a e k e k - l | a.  e k e k k 
k=0  K k=0  k 

where  <J>k  = 0 if  ak  > °»  = * if  < °* 

For  further  reference,  let  us  state  this  result  as 

Lemma  2.2.  If  the  equation  h(p+iv,a,r)  =0  is  satisfied  for 

, -VYk-r 

some  real  u,v,  then  the  lengths  { | a ^ | e , k = 0,1,.. . ,N} 
can  fhrm  a closed  polygon;  that  is,  no  one  of  these  terms  is 
greater  than  the  sum  of  the  others: 
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-UY-T  -UY.-r 

(2.1)  I a j I e 5 l I aJc I e » J - 0,1 N. 

J Wj  K 

Following  Henry  [4],  define  P.  = Pj(a,r)t  j * o,l, 
2.....N,  if  they  exist,  by  the  relations 

-yrr  'pjVr 

(2.2)  I a . | e = l |a J e , j = 0,1 N. 

It  is  easy  to  verify  that  each  p^>pq  always  exist,  are 
unique  and 

(2.3)  PN  = p0  if  N = l*  PN  < P0  if  N - 2’ 

Lemma  2.3.  If  0 < Y^*r  < . . . < y^.r,  then 

Z(a,r)  C [PN(a,  r)  , PQ(a,r)]. 


Lemma  2.2  implies  h(y +i v, a , r)  j-  0 in  either  case,  which 
proves  the  result. 
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The  complete  structure  of  Z(a,r)  is  known  for  the  case 
when  the  components  r are  commensurable.  This  will  be 
stated  as 


Lemma  2.4.  If  ^*1  * ^2  * ' * ' * 


are  commensurable,  that  is, 


rk  = nkB  for  SOme  ^ * 0 and  integers  n^,  k = 1 , . . . ,M , then 

- gX 

h(*,a,r)  is  a polynomial  of  degree  Nn^  in  e , 


Nn 


N 


■ X g 


h (X , a , r)  = a„  I!  (e'^-  rj 

w v=  1 


and 

Z(a,r)  * Z (a , r ) - { ■ g-  ln|rv|,  ^ = 1, 2, . . . ,Nn^} . 
Proof:  Obvious. 

Theorem  2.1.  Z(a,r)  is  continuous  in  a in  the  Hausdorff 
metric.  Also,  if  SC  QR*)N  is  a given  set  and  there  exist 
a < g such  that  Z(a,r)  C (a,g)  for  w e S,  then  there  exist 
a 6 > 0 such  that  Z(b,r)  C (a,g)  for  | b - a | < 6. 


Proof:  From  the  relation 


N -Y  *r 

|hC*,b,r)  - h(X,a,r)|  < l |bk-a.|e  k 

k=0  K K. 


for  any  e > 0,  there  is  a 6 > 0 such  that 


|h(X,b,r)  - hCX, a, r)  | < e for  Re  X c [PN(a,r)  - e,  Pg(a,r)  + e] , 

|b-a|  < 6 

that  is,  h(X,b,r)  - h(X,a,r)  0 as  b -*■  a uniformly  for 

Re  X e tPN(a,r)  - e,  PQ(a,r)  + e]  . 

If  p e Z(b,r)  then  h(X+iv,b,r)  = 0 for  some  v = v(b) . 

If,  in  addition,  b -*■  a,  then  every  limit  point  PQ  of  the  set 

Z(b,rj  as  b -*■  a satisfies  Pq  e Z(a,w)  from  Lemma  .1. 

This  shows  that  6 (Z  (b , r) , Z (a , r) ) 0 as  b -*•  a.  Conversely, 

if  p e Z (a , r ) , then  there  is  a 5 = 5(a)  such  that 
h{P+i5(a)  ,a,r)  = 0.  Therefore,  h(.P+  i£ (a)  ,b ,r ) -*-0  as  b -*•  a 
and  Lemma  2.1  implies  P e Z(b,r).  Thus  6CZ(a,r)  ,Z(b,r))  0 

as  b -*•  a and  the  continuity  of  Z(a,r)  is  proved. 

The  last  statement  of  the  theorem  is  also  a consequence 
of  an  argument  similar  to  the  above. 

Our  next  objective  is  to  discuss  the  dependence  of  ZCa,r) 
on  r.  The  following  example  given  by  Silkowski  [7],  shows 
this  problem  is  much  more  difficult. 

Example  2.1.  Let 

h (X  , r ) = hCX,rlfr2)  = 1 

For  r = Cl, 2),  that  is, 

h(X, 1,2)  = 1 + \ 


•T~  . : > 


-TT  7TCBQ 
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it  is  easy  to  see  that  the  zeros  of  h(X,l,2)  satisfy 
Re  X = -(In  2)/2  < 0.  Therefore,  Z(r)  = {-(In  2 ) / 2 } if 

r = (1,2). 

Now  let  us  consider  r = (r^.r^)  close  to  (1,2).  In 
particular,  take  r = (1-1/ (4n+3) , 2)  where  n is  any  non- 
negative integer. 

It  is  easy  to  verify  that 

h(i (4n+3) ff/2 , 1 - l/(4n+3),2)  = 0. 

Therefore,  Z(r)  D {0}  if  r = (1-1/ (4n+3) , 2)  and  so,  Z(r) 
is  not  continuous  in  r. 

The  numbers  pq»p2  ^or  example  are  Pq  = 0, 

P2  * -In  2 for  r = (1,2).  Also,  Pg(r)  = 0 for  all  r. 
Therefore,  Z(r)  C {P2Cr),0J  where  P2(r)  * 'ln  2 as 

r.-  (1,2). 

What  is  happening  to  the  zeros  of  h,  in  Example  2.1,  as 
r varies?  By  Rouche's  theorem,  for  any  given  rQ  and  any 
compact  set  K in  C for  which  no  zeros  of  h(X,r^)  lie  on 
3K,  there  is  an  e > 0 such  that  | r- rQ | < e implies 
h(X,r)  has  the  same  number  of  zeros  as  h(X,rQ)  in  K. 
However,  a small  change  in  r does  not  necessarily  give  a 
small  change  in  h uniformly  in  a strip  as  was  the  case 
when  the  coefficients  were  varied  as  in  Theorem  2.1.  The  non- 
compactness of  the  strip  plays  an  essential  role  when  r is 
varied. 

For  the  purpose  of  intuition,  it  is  worthwhile  to  note  the 
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following  fact  about  Example  2.1.  For  r = (1,2),  the  zeros 
of  h belonged  to  a vertical  line  Re  X = -(In  2)/2  and  were 
given  by  X = -(In  2)/2  + i(tan  */7  + 2kfi),  k = 0,  1,  2,...  . 
For  a small  change  in  r,  this  vertical  line  of  zeros  is  moved 
a large  distance.  In  fact,  it  may  include  Re  X = 0 . The 
figure  below  is  instructive 


l 

I 

I 

I 

I 


I 

I 

I 

l 

t 

I 

: t 


W = (1,2) 


w = (1  + e , 2) 


We  shall  see  below  that  it  is  actually  possible  for  the 
real  parts  of  the  zeros  of  h to  fill  an  interval. 
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The  above  example  shows  that  Z(r,a)  is  not  necessarily 
continuous  in  r.  However,  it  is  if  the  components  of  r are  rationally 
independent.  This  is  the  content  of  the  next  theorem  where  we  write 
Z(r)  = Z(r,a)  since  a is  fixed. 


Theorem  2.2.  If  rQ  e QR*)M  is  fixed  and  the  components  of  rQ 
are  rationally  independent,  then  Z(r)  ■+  Z(rg)  in  the 
Hausdorff  metric  as  r -*■  r^. 


Proof : Suppose  P(r)  e Z(r),  hCP(r)  + i°Cr),r)  = 0 for  some 
real  o(r)  . If  r -*■  rQ  we  may  assume  P(r)  -*•  PQ. 

Consider  h (Pn+ i v , r n)  . 


h(Viu’ro)  ■ 


r 'Wr0  'ivYr0  . 

L a,  e e 


? 'Vk'ro  •ioW'Vr  iYt°Wr-vr0> 

= l a,  e e e = 

k=0  K 

N -(P0+io(r))Y  -r  PQYk-  (r-rQ)  iYk- (a(r)r-vr0) 
= l a,  e e e 

k=0  K 


By  Kronecker's  Theorem,  for  any  sequence  r^  -*•  r^ , choose 
{v.  },  v.  -v  oo  as  z -*■  » such  that 


iYi,  • (o(r . ) r . - v rn) 
, k ^ v J J J , e 0 


-*1  as  £ -*•  00 , 


By  the  diagonalization  procedure,  we  can  choose  a subsequence 
{v  ),  v - - as  j such  that 
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‘V^VVVo’ 


■}•  1 as  j 


Thus,  h(pQ+i^j,rj)  -*•  0 as  j -*■  ® and  every  limit  point 


Pq  of  Z(r)  satisfies  P^  £ Z(rQ).  This  shows  that 
6(Z(r)  ,Z  CrQ))  -•>  0 as  r - rQ. 


Conversely,  suppose  P £ Z(rg).  Then  there  exist  a cr 
such  that  h(P+io,r^)  = 0.  We  also  have 


N -(P*io)Y kT0  -CPHo)Y.-Cr-r J 

hCf^icr, r)  = I a,e  K ue  k u 

k=0  K 


0 as  r -*■  rQ. 


Therefore,  6 (Z  (r^)  , Z (r) ) -*•  0 as  r -►  r0.  This  proves  the 


theorem . 


As  an  immediate  consequence  of  Theorem  2.2,  we  have  the 
following  result. 

Corollary  2.1.  2T(r)  is  lower  semicontinuous  in  r;  that  is. 


lim  inf  Z(r)  = T(r) 
r-vru 


3. Characterization  of  Z(r,a).  The  following  characterization 

of  Z"(r)  = Z"(r,a)  was  stated  without  proof  by  Henry  [4]  . 


Theorem  3.1.  If 


and  the  components  of  r are  rationally  independent,  then 
P e ZO)  if  and  only  if  there  is  a 9 such  that 
110,0, r)  = 0. 


Proof:  If  h(p+iv,r)  = o,  then  3 9 = vr  such  that 

HCp,9,t)  = 0. 

Conversely,  if  there  exist  9 = 9 e [0,2r] 

j = 1 , . . . ,M,  such  that 

N -py.  • r i y,  * 0 

a0  * l akc  v e ■ 0. 

U k=l  K 


By  Kroneckcr's  Theorem,  there  exists  a sequence  {vn}f 
such  that 

iYk-(0-v"r) 

e -*■  1 as  n -*■  ». 


Therefore, 
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N ■ P Y.  ■ r -ivnY.-r 
h(P+ i^n ,i)  - an  + I a.e  c 

0 k=l  K 

N - PY.  * r - i Yv®  iYv-(9-vnr) 

= an  + I a,e  c e ♦ 0 

0 k-1  K 


as  n -»  °°.  But  this  implies  that  P £ Z(r). 

Theorem  2 . 2 states  that  ZCr)  is  continuous  at  those 
vectors  r with  rationally  independent  components  and 
Theorem  3.1  gives  a way  for  computing  Z(r)  at  such  vectors  r. 


An  important  consequence  of  Theorem  3.1  is  the  following 


result . 


Corollary  3.1.  The  following  statements  are  equivalent 

(i)  0 e Ifr0)  for  some  r°  with  rationally  independent 

components . 

(ii)  0 e Z’(r)  for  all  r with  rationally  independent 
components . 

Proof:  Since  H(0,6,r)  in  Relation  (.3.1)  is  independent  of  r, 

it  is  clear  from  Theorem  3.1  that  (i)  =>  (ii) • The  other  way  is 
obvious . 

Another  easy  consequence  of  Theorem  3.1  is 


Proof ; If  the  components  of  r are  rationally  independent, 
then  Z( r)  is  characterized  by  the  solutions  of  H(p,9,r)  = 0. 

Since  these  solutions  are  analytic  varieties,  it  is  impossible 

to  have  the  following  property:  there  exists  a P e Z(r) 

{Pj}“=1  C Z(r),  Pj  + P as  j -*■  CPj  + 1*Pj)  H Z(r)  = <p . This 
proves  the  corollary  when  the  components  are  rationally  independent. 

For  any  r e 0R*)M  there  exists  a 3 e 0R*)q  for  some  integer 
q such  that  the  components  of  q are  rationally  independent. 

Apply  the  previous  result  to  3 to  complete  the  proof. 

Another  easy  consequence  of  Theorem  3.1  and  Theorem  2.2  is 

Corollary  3.2.  If 

p(r)  = min{p:  39  e 1RM  with  H(p,0,r)  = 0} 

a(r)  = max{p  39  e IR^  with  H(p,0,r)  = 0} 

(3.2)  x (r)  = max{p(r)  < p < 0:  30  e IR^  with  H(p,0,r)  = 0} 

if  p(r)  < 0 

x+(r)  = min{o(r)  > p > 0:  30  e IRM  with  H(p,0,r)  = 0} 

if  o(r)  > 0. 

Then  p(r),  a(r) , x_(r),  x+(r)  are  continuous  in  r,  are  either 
f _Q,  for  all  r or  = 0 for  all  r,  and 

Z(r)  C [P(r),  x (r) ] u (x+ (r) , o(r) ] . 

Furthermore,  p(r),  x_(r),  x+(r),  o(r)  e Z(r) . if  the  components 
of  r are  rationally  independent.  Of  course,  it  is  understood  that 
the  interval  (p(r) ,x_  (r ) ] (respectively,  [t+ (r) , a(r) ] ) is  not 
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considered  if  o(r)  < 0 (respectively,  P(r)  > 0). 

We  remark  that  one  could  give  a finer  structure  theorem  for 
T(r)  than  Corollary  3.2  by  specifying  a finite  of  number  of 
intervals  which  vary  continuously  with  r and  which  coincide  with 
Z^(r)  when  the  components  of  r are  rationally  independent. 
However,  the  number  of  disjoint  intervals  would  not  be  constant 
(examples  will  be  given  later)  in  r.  On  the  other  hand,  the 
structure  theorem  in  Corollary  3.2  is  independent  of  r.  In  fact, 
for  any  r^  with  rationally  independent  components,  there  is  a 
neighborhood  U(r^)  of  r*^  such  that  only  one  of  the  following 
situations  occur: 


(i) 

T.(r)  = T+(r) 

= 0 

for 

all 

r;  that  is,  Z(r)  contains 

zero  for  all 

r e 

U(r°) 

• 

(ii) 

T_(r)  < 0 < t 

for 

all 

r e U(r°);  that  is,  Z(r) 

contains  elements 

< 0 

and 

> 0 for  all  values  of 

r e U(r°)  . 

(iii) 

x_(r)  = o(r) 

(t  (r)  does  not  exist);  that  is  either 

Z(r)  n [0,») 

= * 

for  all 

r e U(r°)  or  Z(r)  n (O,00) 

for  all  r e 

U(r°) 

• 

(iv) 

*+(r)  = p(r> 

(t  (r)  does  not  exist) ; either 

Z(r)  n (-“,01 

= ♦ 

or 

Z(r) 

O (-°°,0]  = (0)  for  all 

r e U(r°) . 

These  remarks  will  be  related  to  stability  in  a later  section. 


4.  A special  case.  When  the  function  h(X,r;a)  has  the  special 


form 


! v 


L 
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(4.1) 


-Ar. 


h(A,r,a)  = 1 + ^=1aj.e  ^ 


one  can  give  a more  precise  description  of  the  set  Z(r)  = Z(r,a). 

This  corresponds  to  the  case  where  N = M and  y.,  =0  if  j ^ k, 

1 K 

Y j j = 1,  j = 1,2,...,M.  It  is  the  purpose  of  this  section  to 
discuss  the  zeros  of  the  function  h in  Relation  (4.1). 


Theorem  4.1.  Suppose  0 < rj  < . . . < rN  and  define  PQ,...,PN 
If  the  set  {r^,  k = 1,2,...,N}  is  rationally 

independent,  then  P e Z(a,w)  if  and  only  if  { | a ^ | , 

-Prv 

| a ^ | e K, , k = 1,2, ...,N}  can  form  a closed  polygon.  Also, 

[p^(a ,r) , Pq (a ,r) ] is  the  smallest  closed  interval  containing 

Z(a,  r)  and  Z(a,r)  is  a finite  union  of  closed  intervals. 

In  fact,  if  I.  C [ PM(a , r)  , p (a , r) ] , j = 1,2,...,N  - 1 is 

J n u -Pr.  -Prk 

the  set  (it  may  be  empty)  such  that  | a - I e J > l |a.|e 

3 k/j  k 

N-l 

for  P £ Ij,  then  Z(a,r)  = tPN , PQ1  \ U^Ij. 

Proof:  If  N = 1,  the  theorem  is  trivial.  Thus,  assume  N > 2 

and  define  a = 1,  r^  = 0, 


-Pr . 


- f.j(P)  = | a . | e J - l 
3 3 k^j 


-Pr, 


ak>e 


, j - 0,1,2,. ..,N. 


-Pr, 


The  set  ( | 3^ | , |ak|e  , k = 1,2, ...,N)  can  form  a closed 
polygon  if  and  only  if  f^(p)  < 0 for  all  j = 0,1,2,...,N.  The 
function  H in  Theorem  3.1  for  (4.1)  is 
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N -Pr  . i.6  . 
H(P,9,r)  = an  + l a.e  3e  3 
0 j=l  3 


It  is  clear  that  "f^fP)  < 0 for  all  j and  some  P"  is  equivalent 
to  ’’there  exist  a 6 e irN  such  that  H(P,9,r)  = 0".  Thus, 
the  first  part  of  the  lemma  is  proved.  The  second  part  is 
Corollary  3.2.  The  last  part  is  simply  writing  down  explicitly  what 
it  means  to  have  f . (P)  > 0 for  some  j.  This  proves  the  theorem. 


Corollary  4.1.  Suppose  0 < ^ < . . . < rN  and  define  Pg,...,PN 
by  Relation  (2.2).  If 

pct(r)  = max{  Pj  (r)  < 0,  j = 0,1,2, ...,N) 

Pe(,r^  = min{ Pj  (r)  > 0,  j = 0,1,2,.  ...N) 


(one  of  these  numbers  may  not  exist) , then  pa(r)*  P6(r)  ,re  con- 
tinuous  in  r and  are  either  ^ 0 for  all  r or  = 0 for  all  r. 
Furthermore, 

*<r)  £ IPN(r),Pa(r)l  U lP6Cr).P0(r)l 

and  the  end  points  of  these  intervals  belong  to  T(r)  if  the 
components  of  r are  rationally  independent.  Of  course,  it  is 
understood  the  interval  is  not  considered  if  an  endpoint  does  not 
exist . 

Proof:  This  is  a consequence  of  Corollary  3.2. 
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5.  Stability  and  hyperbolicity.  In  this  section  h(A,a,r)  is 
the  function  defined  in  Relation  (1.1);  that  is 


(5.1) 


h(*,r,a)  - 1 * Ej„ia](e 


-Avr 


We  need  the  following  definitions. 

Definition  5.1.  The  function  h(A,r,a)  is  said  to  be  hyperbolic 
at  r°  if  0 f.  Z"(r®,a).  The  function 

h(A,r,a)  is  hyperbolic  locally  at  r^1  if  there  is  a neighborhood 
U of  r^  and  6 > 0 such  that  Z(r,a)  D [-6,6]  = <j>  for  all 

r e U.  The  function  h(A,r,a)  is  hyperbolic  globally  in  r if 

+ M 

0 i-  Z(r,a)  for  each  r e QR  ) . 


Definition  5.2.  The  function  h(A,r,a)  is  said  to  be  uniformly 
asymptotically  stable  at  r°  if  h(A,r^,a)  is  hyperbolic  and 
2f(r^,a)  n tO,00)  = <t> . It  is  uniformly  asymptotically  stable  locally 
at  r^,  if  it  is  hyperbolic  locally  at  r^  and  Z’(r.a)  H ( 0 ,°°)  = <t> 
for  r e U.  It  is  uniformly  asymptotically  stable  globally  in  r 
if  it  is  hyperbolic  globally  in  r and  T(r,a)  H CO,00)  = $ 
for  all  r. 

We  now  prove  the  following  fundamental  result.  In  the  statement 
of  the  theorem,  P(r),o(r),x  (r),t+(r)  are  defined  in  Relation  (3.2) 


r 
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Theorem  5.1.  The  following  statements  are  equivalent. 

(i)  There  is  an  re  QR+)M,  r = (r^...,^),  with  the  set 
{ r ^ rationally  independent,  such  that  the  function  h(A,r,a) 

is  hyperbolic  at  r . 

(ii)  h(A,r,a)  is  hyperbolic  locally  at  some  r®. 

(iii)  h(A,r,a)  is  hyperbolic  globally  in  r. 

(iv)  There  is  an  r°  e QR+)^  and  a neighborhood  U of  r9 

such  that  h(A,r,a)  is  hyperbolic  for  every  r e U with 
the  components  of  r commensurable. 

(v)  t (r°)  < 0,  T+(r°)  > 0 for  some  r°  e QR+)M  if  these 
numbers  exist. 

(vi)  If 

m -Ar . 

(5.2)  h(A ,r,a)  = det[I  - F?=1A.e  3 ] 

then 

\ 

m - i0  • 

C5.3)  (w(0):  dettPl  - Ij=1Aje  3]  = 0,  6 e]RM}n{|y|  = 1}  = 4>. 

Proof : Let  us  first  prove  (v)  <=>  (vi) . If  h(A,r,a)  is  given  by 
Relation  (5.2),  then  the  function  H(A,0,r)  in  Relation  (3.1)  is 
given  by 

M *Pri  i9-i 

H(P,0,r)  = det  [I  - Ij=1Aje  3 e 3 ] . 

The  equivalence  of  statements  (v)  and  (vi)  is  now  immediate. 

From  the  definitions  of  * (r),x+(r)  and  the  remarks  following 
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Corollary  3.2,  we  have  (v)  < = > (i) , (v)  < = > (iii)  < = > (iv)  . 

Obviously  (iii)  <=>  (iv) . To  complete  the  proof  of  the  theorem, 
we  show  (iv)  <=>  (v) . If  (iv)  is  satisfied  and  (v)  is  not,  then 
t (r)  = i+(r)  = 0 for  all  r from  the  remarks  after  Corollary  3.2. 
Since  T_(r),i+(r)  e Z”(r)  if  the  components  of  r are  rationally 
independent  and  Z^r)  is  continuous  in  r at  these  points,  this 
is  an  obvious  contradiction.  This  proves  the  theorem. 

Since  stability  is  so  important  in  the  applications,  we  re- 
state Theorem  4.1  for  this  case. 

Theorem  5.2.  The  following  statements  are  equivalent 

(i)  There  is  an  r e 0R+)M,  r = (rj,...,rJ  with  the  set 
M 

{rj}j_j  rationally  independent,  such  that  the  function  h(A,a,r) 
is  uniformly  asymptotically  stable. 

(ii)  h(A,a,Tp)  is  uniformly  asymptotically  stable  locally  at 
some  rQ  e QR  ) . 

(iii)  h(A,a,r)  is  uniformly  asymptotically  stable  globally 
in  r . 

(iv)  There  is  an  r^  e (IR+)M  and  a neighborhood  U of  r^ 
such  that  h(A,a,r)  is  uniformly  asymptotically  stable  for  every 
r e U with  the  components  of  r commensurable. 

(v)  c(r)  < 0 for  some  r e (R+)^. 

(vi)  If 


h(A , r , a) 


det  [I 


yM  a 
zj=l  j 


-Ar. 

J 


] 


then 
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sup{|u(9)|:  det[y(0)I  - ^=1A.e~  J],  0 e]RM}  < 1. 


i0 


>Mi 


Historically,  Theorem  5.2  developed  in  the  following  way. 
Melvin  [6]  proved  the  result  for  the  scalar  equation  where 
Condition  (vi)  becomes  the  simple  condition 


-N 


Ij.ilAjl  < 1.  A.  eK. 


Hale  [3]  proved  (ii)  <=>  (iii)  in  the  general  case.  Silkowski  [7] 
introduced  the  equivalent  conditions  (i)  aid  (vi) . 


6.  Examples.  In  this  section  we  collect  some  examples  to 
illustrate  the  above  results.  Throughout  the  section,  the  numbers 


P . (r)  are  defined  in  Relations  (2.2),  the  numbers  P(r),  o(r)  , 


T_(r),T+(r)  in  Relations  (3.2) 


Example  6.1.  Let  us  reconsider  Example  2.1;  that  is,  the  function 


1 i -Xr2 

h(X,r)  = 1 + j e 1 + | e 


We  have  seen  that  Pq(t)  = 0 f°r  aH  r an^»  for 


rQ  ■ (1>2),  Z(rQ)  = {-(In  2 ) / 2 } . Now,  Theorem  4.1  implies 
that,  for  any  r=  (r^.r^)  with  rj>r2  rationally 
independent,  [p2  (rl  , Pp(r ) ] s [ P2 Cr) » 0 ] is  the  smallest 
closed  interval  containing  Z(r)  and  P2(r)  is  continuous 
in  r,  p 2 C rQ ) = -In  2.  Furthermore, 


f 


< II 


. 


i 


! 


* 
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1 -P.r2  * Pr , 

2 Ce  - e ) + 1 > 0,  for  P < 0,  ri  < r2 

Therefore  Lj  of  Theorem  4.1  is  the  empty  set  and  Z(r)  = 

[ P 2 Cr ) , 0 ] . Thus,  for  any  neighborhood  U of  rQ,  there  is 
an  r e U such  that  Z(r)  is  a complete  interval  of  length 
approximately  In  2 whereas  for  r = r^,  Z(.rQ)  is  a single 
point . 

Example  6.2.  As  for  Example  6.1,  one  shows  that 
Z(r^)  = l P 2 » Pq  ^ ~ I'*27,.37J  for  the  function 

h(A,r°)  = 1 + e'X  + e'77*,  r°  = (l,n). 

Example  6.3.  Consider  the  equation 

-Ar1  -Ar, 

(6.1)  h (A  , r , a)  = 1 + a^  + a2e  = 0 

where  0 < rj  < r2  and  aj>a2  are  reai  constants.  The  numbers 
Pj , j = 0,1,2,  are  defined  by 


P0r2 


I a2 I e = 1 ' 1 ai 1 e 


'P0rl 


I a2  I e 


I a?  I e 


Plr2 


P2r2 


"plr  1 

I a! I e - 1 


1 + | a e 


‘p2r2 


(6.2) 
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As  remarked  earlier,  < Pq.  The  constant  P^  may  or  may 
not  exist.  From  (6.2),  it  is  clear  that  PQ  < 0 iff 
la^l  + | a 2 1 < 1.  Thus,  h(X,r,a)  is  uniformly  asymptotically  stable 
globally  in  r if  and  only  if  | a 2 1 > 1 + |a^|.  This  means 
h(X,r,a)  is  hyperbolic  globally  in  r and  has  Z(r)  n (-°°,0]  = 4> 
if  and  only  if  | a 2 1 > 1 + |a^|. 

Let  us  now  analyze  the  other  regions  in  the  (a^,a2)  parameter 
space.  If  | a 2 | < 1 + | a^ | (that  is,  P 2 (r ) < 0)  and 
| a^ | + | a2 | > 1 (that  is,  PQ (x)  > 0),  then  we  know  from  the 
general  theory  that  t (r)  < 0 (which  automatically  implies 
t+(r)  > 0)  is  equivalent  to  hyperbolic.  The  number  x_(r)  can 
be  related  to  the  solutions  P^  of  liquation  (6.2).  In  fact,  from 
Theorem  4.1,  x _ (r ) <0  if  and  only  if  there  is  a solution 
Pll(r)  of  Equation  (6.2)  satisfying  p2(r)  < pu(r)  < 0 and 

■plr?  'piri 

| a2 1 e < | a jJ  e - 1 for  Pj^r)  < px  < 0. 

If  | a 2 1 < | a ^ | - 1,  then  there  is  a P^^r)  satisfying  the  above 
properties  and  h(X,r,a)  is  hyperbolic  with  Z(r)  O (-®,0)  f , 

Z(r)  O (0,®)  t <P. 

If  | a 2 1 > |a1|  - 1 


;pir2 


> 


‘plrl 


1 for  Ip^I  < 6 


then 
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for  some  6 > 0,  If,  in  addition,  |a^|  + | a ^ I > 1,  | a £ I < 1 + |a^| 
then  0 e Z”(r)  by  Theorem  4.1  if  the  components  of  r are 
rationally  independent.  Thus,  the  function  h(X,r,a)  is  not 
hyperbolic  globally  in  r. 

In  summary, 

(i)  h(X,r,a)  is  uniformly  asymptotically  stable  globally 
in  r if  and  only  if  |aj|  + |a2|  < 1. 

(ii)  h(X,r,a)  is  hyperbolic  globally  in  r with 

Z(r)  D (-°°,0)  ? <P  if  and  only  if  |a2|  > 1 + | a1 1 . 

(iii)  h(X,r,a)  is  hyperbolic  globally  in  r with 

T(r)  D (-»,0)  f 4> , Z"(r)  n (0,°°)  ? <P  if  and  only  if 

I ai I > 1 + la2l  • 

(iv)  h(X,r,a)  is  not  hyperbolic  globally  in  r if  the  co- 
efficients ai>a2  do  not  satis£y  one  o£  conditions 
in  (i)-(iii). 

The  structure  of  the  set  Z(r,a)  obviously  changes  as  the 
parameter  a varies  from  the  region  in  case  (iii)  to  the  region 
in  case  (ii)  above  since  two  intervals  had  to  merge  as  X(r,a) 
moved  to  positive  axis.  This  structure  can  also  change  even  when 
the  parameters  always  stay  in  a region  corresponding  to  one  case. 
In  fact,  suppose  |a^|  + |a2|  < 1;  that  is,  uniform  asymptotic 
stability  globally  in  r.  Since  | a^ | - 1 < 0,  there  is  an  a2 
sufficiently  small  so  that  the  equation 


-Pr-  -prl 

I a 2 1 6 =laile  1 


has  two  distinct  negative  solutions  P^Cr)  < Pl2(r) 


25 
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[ P 2 (r) > Pq(t)] • Theorem  4.1  allows  one  to  conclude  that  T(r,a) 
consists  of  two  intervals. 

Let  us  make  one  other  remark  about  this  example.  The  number 
of  intervals  in  Z(r,a)  may  also  change  with  r.  In  fact, 
suppose  | a 2 1 = | | - 1.  The  function  h(A,r,a)  is  not  hyper- 
bolic globally  in  r in  this  case.  The  equation 

Ant  -Pr0  -pr. 

f(P,a,r)  d=f  i a 2 | e - laje  +1  = 0 

has  the  solution  P = 0.  Since 


9f C0,a,r)  _ 
5p 


= -\a2\r2  + \a1\r1 


and  f(P,a,r)  -*-1  as  P -*•  °°s  there  will  be  a positive  zero 
of  f if  I a 2 1 r 2 51  lai_lrj'  Since  f(P,a,rl  -+  +0°  as 
p ■+  *"  there  will  be  a negative  zero  of  f if  ^2^2  < 

I a 1 1 • 

Therefore,  if  l a 2 ^ r 2 ^ 1 a 1 lr 1 * that  *s» 
r2  t C|a2|/(l+|a2|))  r2»  ri»r2  rationally  independent,  the 
set  T(a, r)  will  consist  of  two  intervals.  When 
ri  “ C | a2  J / Cl-*-  |a2  D)  r2  the  set  Z(a,r)  will  consist  of  one 
interval. 

T 

j 

t 

Example  6 . 4.  Consider  the  equation 


f 


Let  us  study  ZCe)  as  e 0 and  always  assume  that 
C2  > C1  > °* 

As  a first  case,  if  e2  = 2 £j  then  h(X,E)  = 0 if  and 
'XE1  1 

only  xf  1 - ce  = 0,  Re  A = A-  ln|c|.  Thus,  if  |c|  > 1, 

Re  X +®  as  -*•  0;  if  |c|  = 1 , Re  X = 0 for  all  e^; 
if  |c|  < 1,  then  Re  X + -«  as  e^  -*•.(). 

If  e2  51  E1  > °’  we  know  that  z(c)  c l p2 ^ > pq ( £) J where 

P2  = P2^C^’  p0  = P0^>  satisfy  the  equations 


(6.4) 


(a)  1 - 2lc|e’eiP°  ♦ cV2'0 


-e.p. 


7 7 y 7 

Cb)  c e 1 = 1 + 2| 


-c,  P- 


Now  suppose 


(6.5) 


-|c|  < < | c | . 


If  relation  C6.^  ) is  satisfied,  then  2 | c | + c >1  and 
p0  = p CO  > 0,  p0(E)  -*■  +«  as  e -*■  0.  Furthermore,  if 
P2  > 0,  then 


■^•jPp  J ~^7^7  7 _£1P? 

1 ♦ 2 | c | e 1 2 = c2e  2 2 < c2e  12 

7 _SP  7 

1 < (c  -2  | c | )e  1 < c2  - 2 | c | 


< - c . 
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Thus,  if  Relation  (6.5)  is  satisfied,  then  ^ = e2^  < 

P2Ce)  -•  as  e 0. 

Also,  if  e2  > £l  > 0 are  rationally  independent,  then 
(P2Ce) » PQ(e) 1 is  the  smallest  interval  containing  Z(£). 
Thus,  if  Relation  (6.5)  holds,  the  smallest  closed  interval 
containing  Z(£)  approaches  C-®,*00)  as  £ -*■  0. 

To  determine  when  Z(£)  is  a single  interval,  we  should 
find  P1(£).  The  number  P^Ce),  if  it  exists,  must  be  a 
zero  of  the  function 

? •pe?  'pei 

fCP,c)  = c e 2 - 2 | c | e 1 * 1 


this  function  has  a unique  minimum  at  a point  given  by 


a 


In 


If  |c|  <1,  then  we  can  choose  ^ > such  that 
| c | ^2^ 2 el  = * and  tflUS  a “ 0*  Since  f(0,£)  = (|c|-l)^>  0 
if  |c|  < 1,  it  follows  that  f(P,£)  > 0 for  all  p and  P^ 
does  not  exist.  This  means  that  Z(a,£)  = lP2 (e) » Pq (e)  1 • 

We  can  thus  choose  ej»e2  ei  < e2’  SO 

\z\z2/2t-1  = 1 and  Z(a,£)  -*■  (*®,«). 

If 


then  P0(E)  < 0 and  Pg(E)  *" 


as  c ■*  .0.  If 
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(6.7) 


Cl  * 


then  PQ(e)  = 0 for  all  e2  > ei  > °»  p2^e^  < °*  p2^e^ 
as  e -*•  0 and  the  smallest  closed  interval  [P2(c),0], 
containing  Z(e)  approaches  (-®,0]  as  e 0. 


Example  6.5.  I.ct  r]<r2<rj»hC^,r)  = l + c + 

- X r - Xr  .. 

e L + c p.  If  (r1,r2>r.)  are  rationally  independent, 
then  the  smallest  closed  interval  containing  Z(r)  is 
(P3,P0)  where 


•P.r.  *P-ri  * P7 

3 3 . a 1 , 3 

e = 1 + e + c 


i = c 


'Vi 


'P0r2  _P0r3 

+ e + c 


The  numbers  » p2  arc  defined  by 


;Plrl  . x + e"Plr2  + e‘Plr3 
;P2r2  . ! + e*  P2rl  + c-p2r3 


if  they  exist.  This  implies  necessarily  that  Pj  <0,  P2  < 0. 
On  the  other  hand,  for  p < 0,  the  functions 


- Pr^  -pr  -pr 

f(p)  = e J ♦ 1 - e 1 ♦ e L 

-Pr,  -Pr.  - pr 

g(P)  = e ♦ 1 ♦ e 1 - e L 
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are  decreasing  and  positive  for  p = 0.  Thus,  ^o  n°t 

exist  and  Z(r)  = [Pj.PqI  if  (rl,r2,r3^  arc  rationally 
independent . 

Suppose  (rj,r2,r3)  arc  no  longer  rationally  independent; 
in  particular,  r^  = 1 , r^  = 2,  r^  = r.  Then  [pj,Pg]  = 

[-.56,  .60].  IVhat  is  the  smallest  interval  ^ p5  * po  ^ 
containing  T(r0),  rQ  = (1,2,*)? 

Front  Theorem  3.1,  we  need  only  determine  0 = anc*  p 

such  that 


n i9,  ,n  2i9.  _o  i9 
H(p,9)  = 1 + e"Pe  1 + e‘2pe  1 + c pe  2 = 0 


that  is, 


Re  X 


These  curves  intersect  if  ancl  only  if  P e [p,,o]  where 


Pj  is  as  above  and  cr  satisfies 


*"o  /j  ( , -2c. 

c - -T  Cl-e  ). 


Thus,  Z(rQ)  = I P ^ » 0 1 * I-.56, .30)  and  [ , a]  CjlPjiPgJ* 

There  is  a definite  shrinking  of  the  interval  when  the  delays 
are  not  allowed  to  vary  independently. 
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